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ABSTRACT 


This  thesis  is  part  of  a  continuing  study  in  developing  convective  heat  transfer 
correlations  for  a  cylinder  in  a  high  amplitude  zero-mean  oscillating  flow.  The 
experiment  described  here  utilizes  the  RTD  technique  and  a  steady  state  heat  transfer 
measurement  method  with  a  platinum  wire,  serving  as  the  test  section,  positioned  across 
the  inner  diameter  of  a  cylindrical  plexiglass  chamber  supporting  a  strong  resonant  axial 
acoustic  field.  Utilizing  two  different  vsire  diameters  of  0.050  mm  and  0.127  mm, 
various  pressure  ratios,  frequencies,  and  temperature  differences,  separated  flow  heat 
transfer  correlations  have  been  developed.  This  work  would  find  application  in  the 
design  of  heat  exchangers  for  a  thermoacoustic  engine. 


V 


VI 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

II.  HISTORICAL . 3 

A.  REFRIGERATION . 3 

1 .  Vapor  Refrigeration  Cycle . 3 

2.  Thermoacoustic  Cycle . 4 

B.  HEAT  TRANSFER . 6 

1 .  Natural/Forced  Convection . 6 

2.  Thermoacoustic  Streaming . 8 

in.  FLUID  MECHANICS . 11 

A.  RADIAN  WAVELENGTH . 11 

B.  DISPLACEMENT  AMPLITUDE . 12 

C.  FREQUENCY  PARAMETER . 14 

D.  REYNOLDS  NUMBER . 15 

E.  FLOW  SEPARATION . 16 

IV.  EXPERIMENT . 19 

A.  PURPOSE . 19 

B.  EQUffMENT . 20 

1.  Test  Cylinder . 20 

2.  Acoustic  Chamber . 21 

3.  Electronics . 22 

C.  PROCEDURE . 25 

V.  RESULTS  AND  DISCUSSION . 31 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS . 39 

APPENDIX  A:  CALIBRATION . 43 

APPENDIX  B:  SAMPLE  CALCULATION . 49 

A.  NUSSELT  NUMBER . 49 

vii 


B.  REYNOLDS  NUMBER . 51 

APPENDK  C:  UNCERTAINTY  ANALYSIS . 53 

A.  PLATINUM  WIRE  RESISTANCE . 53 

B.  NUSSELT  NUMBER . 55 

C.  REYNOLDS  NUMBER . 57 

APPENDDCD:  EXPERIMENTAL  DATA . 61 

LIST  OF  REFERENCES . 95 

INITIAL  DISTRIBUTION  LIST . 97 


viii 


LIST  OF  FIGURES 


Figure  1.  Basic  Vapor  Refrigeration  Cycle . 3 

Figure  2.  Basic  Thermoacoustic  Refrigeration  Components . 5 

Figure  3.  Vertical  Cylinder  Natural  Convection . 7 

Figure  4.  Acoustic  Streaming  Pattern . 9 

I 

Figure  5.  Isothermal  “Streaked”  Flow . 17 

Figure  6.  Acoustic  Chamber  and  Accessories . 22 

Figure  7.  Acoustic  Electronics  Package . 23 

Figure  8.  Heat  Transfer  Electronics  Package . 24 

Figure  9.  Standing  Wave  Velocity  Profile  in  a  Circular  Cylinder . 26 

Figure  10.  Nusselt  Number  versus  Reynolds  Number  (d  =  50.8  pm) . 35 

Figure  11.  Nusselt  Number  versus  Reynolds  Number  (d  =  127  pm) . 36 

Figure  12.  Nusselt  Number  versus  Reynolds  Number  (both  wires) . ..37 

Figure  13.  Curve  Fit  for  Nusselt  Number  versus  Reynolds  Number  (d  =  50  pm) . 41 

Figure  14.  Calibration  Chamber . 44 

Figure  15.  Platinum  Wire  Calibration  Curve  (d  =  50.8  pm) . 46 

Figure  16.  Platinum  Wire  Calibration  Curve  (d  =  127  pm) . 47 


IX 


X 


LIST  OF  TABLES 


Table  1.  Platinum  Wire  Calibration  Data  (d  =  50.8  |a.m) . 45 

Table  2.  Platinum  Wire  Calibration  Data  (d  =  127  |j.m) . 45 


XI 


LIST  OF  SYMBOLS,  ACRONYMS,  AND/OR  ABBREVIATIONS 


a  test  cylinder  radius  [m] 

A  particle  displacement  [m] 

A  particle  displacement  amplitude  [m] 

As  surface  area  [m^] 

c  speed  of  smmd  [m/s] 

d  diameter  of  test  cylinder  [m] 

f  frequency  [Hz] 

h  convective  heat  transfer  coefficient  [W/m^-K] 

HX  heat  exchanger 

I  current  [Amps] 

k  thermal  conductivity  [W/m-K] 

KC  Keulegan-Carpenter  number 

1  length  of  test  cylinder  [m] 

L  distance  from  test  cylinder  to  chamber  end  plate  [m] 
M  Mach  number 

Nu  Nusselt  number 

Pm  mean  ambient  pressure  [Pa] 

Pref  reference  pressure  [Pa] 

PR  pressure  ratio 

R  gas  constant  [J/kg-K] 

Re  Reynolds  number 

RpR  precision  resistor  resistance  [ohms] 

Rs  streaming  Reynolds  number 

Rw  wire  resistance  [ohms] 

RMS  root  mean  square 

RTD  resistance  temperature  detector 

SPL  soimd  pressure  level 

t  time  [s] 

Tc  cold  sink  temperature  [K] 

Th  hot  sink  temperature  [K] 

Ts  surface  temperature  [K] 

Too  ambient  temperature  [K] 

U  velocity  [m/s] 

Uo  velocity  amplitude  [m/s] 

Vm  pressure  transducer  voltage  [volts] 

Vo  voltage  [volts] 

VpR  precision  resistor  voltage  [volts] 

Vw  platinum  wire  voltage  [volts] 

(3  amplitude  ratio 


xiii 


X  cylinder  length  scale 

6  boundary  layer  thickness  [m] 

s  amplitude  parameter 

7  ratio  of  specific  heats 

X  wavelength  [m] 

Z  radian  wavelength  [m/rad] 

A  frequency  parameter 

V  kinematic  viscosity  [m^/s] 

(D  radian  frequency  [rad/s] 


XIV 


I.  INTRODUCTION 


Pollution  free  refrigeration  is  a  concept  worthy  of  close  inspection  in  this  day  and 
age  of  an  environmentally  conscious  public.  A  more  far-reaching  concern  to  the  masses 
than  the  purported  global  environmental  awareness  may  be  the  economic  impact  from  a 
recent  international  agreement  to  ban  all  consumption  of  chloroflourocarbons  (CFCs)  by 
the  turn  of  the  century.  This  will  affect  everyone  in  the  industrialized  world. 
Technologies  are  being  investigated  to  supplant  or  improve  present  day  refrigerants 
and/or  their  associated  cycles.  Researchers,  whose  interests  come  from  a  broad  spectrum 
including  academia,  space  and  the  automotive  industry,  seek  to  provide  an  alternative  to 
today’s  refrigeration  technology  whereby  a  CFC-free  environment  remains. 

One  method,  thermoacoustic  refrigeration,  blends,  in  part,  the  disciplines  of 
convective  heat  transfer  and  acoustics.  Four  real  world  applications  of  this  technology 
include  the  Space  ThermoAcoustic  Refrigerator  (STAR)  [Ref  1],  a  food  refrigerator  built 
in  the  Republic  of  South  Africa,  a  cryogenic  refrigerator  for  electronics  cooling  and  a 
device  built  by  the  Ford  Motor  Company  for  potential  automobile  uses.  Through 
continued  innovation  and  experimentation  these  highly  specialized  designs  may  soon 
become  distant  relatives  of  a  new  age  of  environmentally  safe,  mass  produced, 
commercial  thermoacoustic  refrigeration  units. 

If  theory  has  been  put  into  practice,  why  not  close  the  book  on  laboratory 
experimentation  in  this  field?  At  present,  thermoacoustic  refrigerator  technology  depends 
on  efficient  solid  to  gas  heat  transfer  from  its  heat  exchangers.  Lee  and  Richardson  [Ref 
2]  reiterate  what  has  been  known  for  a  long  time,  that  is,  heat  transfer  coefficients 
between  gases  and  solid  surfaces  are  small.  If  the  turbulence  intensity  of  the  gas  free 
stream  is  raised,  heat  transfer  rates  will  increase.  Researchers  have  found  that 
experiments  and  analysis  involving  stationary  sound  fields  rather  than  real  turbulent  flow 
are  much  simpler.  Their  experiments  have  shown  that  under  favorable  circumstances  the 
presence  of  oscillations  can  cause  significant  increases  in  gas  heat  transfer  coefficients. 
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Today’s  thermoacoustic  refrigerators  lack  the  efficiency  of  modem  day  vapor 
compression  cycles.  Operating  at  best  at  about  50%  efficiency  of  conventional 
refrigerators,  there  is  much  room  for  improvement  especially  in  heat  exchanger  design. 
The  intent  here  is  to  extend  the  experimental  basis  for  heat  exchanger  design  in 
thermoacoustic  applications.  The  fundamental  flow  regimes  and  the  results  from  varying 
characteristic  parameters  within  each  have  not  yet  been  completely  analyzed.  The 
motivation  for  this  work  is  to  develop  a  better  understanding  of  the  relationship  between 
high  amplitude  resonant  sound  fields  and  convective  heat  transfer  rates  in  a  gaseous 
medium  and  to  develop  useful  correlations  for  a  portion  of  the  wide  parameter  regime  of 
interest. 
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II.  HISTORICAL 


A.  REFRIGERATION 

1.  Vapor  Refrigeration  Cycle 

Vapor  compression  systems  are  the  most  common  refrigeration  cycles  used  for 
both  industrial  and  commercial  refrigeration  (Figure  1).  The  working  fluid,  or  refrigerant, 
in  these  systems  is  typically  a  halogenated  hydrocarbon.  A  measure  of  efficiency  for  these 
cycles,  designated  the  coefficient  of  performance  (COP),  comes  fi’om  the  ratio  of  actual 
refrigerating  effect  to  the  net  work  required.  In  theory,  the  Carnot  cycle  represents  the 
most  efficient  cycle.  However,  several  impracticalities  arise,  therefore,  actual  systems 
sacrifice  efficiency  for  a  reduction  in  costs  and  maintenance. 


Figure  1.  Basic  Vapor  Refrigeration  Cycle 


Today’s  refrigerators  grew  out  of  years  of  theoretical  and  experimental  studies  to 
become  efficient,  reliable  and  versatile  machines  serving  numerous  applications. 
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Operating  efficiencies  on  the  order  of  over  90%  of  Carnot  efficiency  are  not  uncommon. 
Recent  international  laws  requiring  an  end  to  the  production  of  environmentally 
hazardous,  ozone-depleting,  hydrocarbon  refrigerants  (i.e.,  R-11,  R-12,  R-114)  used  in 
many  of  these  systems  led  to  new  temporarily  acceptable  substitutes  (i.e.,  R-134a)  and  an 
interest,  at  least  academically,  in  radically  new  methods  for  a  future  standardized  means 
of  refrigeration. 

2.  Thermoacoustic  Cycle 

Thermoacoustic  refrigeration  systems  can  be  described  in  the  following  very 
simplistic  way.  There  are  six  basic  system  components  illustrated  in  Figure  2.  They 
include  the  driver  or  loudspeaker,  two  heat  exchangers,  a  “stack,”  and  a  sound  medium 
(typically  a  gas)  all  enclosed  within  an  acoustic  chamber.  A  resonating  sound  wave, 
contained  within  the  acoustic  chamber,  establishes  a  sinusoidal  varying  pressure  field 
with  antinodes  at  the  chamber  ends.  Acoustic  theory  explains  the  simultaneous  formation 
of  sinusoidal  displacement  and  velocity  fields  90  degrees  out  of  phase  with  the  pressure 
field.  They  are  anchored  by  nodes  at  the  chamber  ends.  The  sound  wave  remains 
stationary,  however,  the  gas  molecules  within  the  chamber  rapidly  move  back  and  forth. 
By  precisely  positioning  the  heat  exchange  elements  (i.e.,  hot  heat  exchanger,  and  cold 
heat  exchanger  with  a  “stack”  of  low  thermal  conductivity  plates  in  between)  away  from 
the  pressure  and  displacement  nodes  a  “bucket  brigade”  means  of  heat  transfer  between 
gas  particles  is  created  along  the  length  of  each  plate  within  the  stack.  If  we  assume  that 
a  half-wavelength  sound  wave  exists  within  the  chamber  of  Figure  2,  the  stack  cools  at 
the  end  facing  the  driver  and  heats  up  at  the  other  end.  The  magnitude  of  the  temperature 
gradient  is  directly  related  to  the  pressure  ratio  created  in  the  chamber.  As  the  pressure 
amplitude  increases  so  does  the  temperature  difference. 

Figure  2  also  illustrates  a  simplified  version  for  one  cycle  of  particle  oscillation, 
[Swift  Ref.  3].  A  gas  particle  undergoes  adiabatic  compression  when  displaced  to 
position  A.  Following  the  laws  of  thermodynamics,  an  accompanying  temperature  rise 
leaves  the  particle  hotter  than  the  plate.  Heat  is  rejected  to  the  plate,  from  the  particle. 
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which  returns  to  position  B  while  undergoing  adiabatic  expansion.  A  temperature 
decrease  accompanying  particle  expansion  lowers  its  temperature  below  that  of  the  plate. 
Heat  flows  to  the  particle  and  “up”  the  temperature  gradient  with  subsequent  compression 
and  displacement  of  the  particle  back  to  position  A.  A  shallow  positive  temperature 
gradient  forms  across  the  stack  from  the  cold  heat  exchanger  to  the  hot  heat  exchanger 
through  adjacent  particle  heat  transfer  across  the  plates  as  described  above. 


Figure  2.  Basic  Thermoacoustic  Refrigeration  Components 

In  1993  A.  D.  Little,  Inc.,  performed  a  research  needs  assessment  for  the 
Department  of  Energy  [Ref.  4].  They  analyzed  more  than  fifteen  refrigeration  alternatives 
including  substitute  refngerants,  improved  components  and/or  alternate  cycles.  For 
domestic  refrigeration,  Little’s  report  concluded  that  thermoacoustic  refrigeration  was  an 
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unlikely  replacement  to  vapor  compression.  Under  the  categories  of  priority  for  research 
and  development,  probability  of  success,  and  industry  involvement,  thermoacoustics 
ranked  lowest  compared  to  all  other  technologies  considered.  Since  that  report,  however, 
several  corporate  sponsored  projects,  mentioned  above  in  part,  affirm  the  viability  of 
thermoacoustic  refrigeration.  Steven  L.  Garrett  et  al.  [Ref.  5]  succinctly  outline  the 
benefits  and  predominant  technical  issues  surrounding  this  technology. 

A  simple  description  of  the  basic  flow  of  heat  from  one  end  of  the  stack  to  the 
other,  in  the  presence  of  an  acoustic  standing  wave,  is  given  above.  What  has  not  been 
considered  in  great  detail  is  the  effect  that  the  same  sound  wave  will  have  on  the  heat 
transfer  from  the  heat  exchanger  to  or  from  the  stack  ends.  If  a  gap  in  physical  contact 
exists  between  a  heat  exchanger  and  the  end  of  the  stack,  how  can  we  best  operate  and 
design  the  refiigerator  to  maximize  efficient  heat  transfer  across  that  gap?  It  is  here  that 
significant  advances  in  this  technology  need  to  be  made  both  in  understanding  the  theory 
and  incorporating  experimental  results  into  design.  This  study  attempts  to  close  the  gap 
between  the  theory  behind  thermoacoustic  refngerator  heat  exchange  in  a  soxmd  field  and 
an  efficient  heat  exchanger  design  by  providing  heat  transfer  correlations  in  the  separated 
flow  regime.  It  is  here  that  great  advances  in  system  efficiency  and  performance  can  and 
need  to  be  made  in  order  to  persuade  hesitant  academics  and  conservative,  potential 
corporate  sponsors  of  the  need  for  future  research,  development  and  investment  in  this 
relatively  new  technology. 

B.  HEAT  TRANSFER 

1.  Natural/Forced  Convection 

Generally  speaking,  there  are  two  distinct  modes  of  convective  heat  transfer, 
natural  (free)  and  forced.  Each  represents  a  different  means  by  which  fluid  motion  is 
sustained.  The  importance  of  one  mode  relative  to  the  other  is  governed  by  the  strength 
of  its  respective  fluid  motion.  Forced  convection  fluid  motion,  or  fluid  velocity,  is 
generated  by  some  external  forcing  condition  and,  hence,  the  flow  and  thermal  fields  are 
largely  imcoupled.  As  an  example,  consider  a  horizontal  cylinder  immersed  in  air. 
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Assuming  that  a  temperature  difference  exists  between  the  cylinder  and  air,  a  fan  used  to 
force  air  over  the  cylinder  creates  the  necessary  fluid  motion  for  forced  convection  heat 
transfer.  On  the  other  hand,  natural  convection  fluid  motion  arises  when  a  body  force 
acts  upon  a  density  difference.  This  density  gradient  typically  results  from  a  temperature 
difference  between  the  body  and  the  medium.  The  result  is  coupled  flow  and  thermal 
fields,  a  more  difficult  problem  to  solve.  Using  a  vertical  cylinder  in  a  quiescent  air 
medium,  a  temperature  difference  between  the  cylinder  and  air  gives  way  to  buoyant 
forces  that  create  natural  convection  currents  adjacent  to  the  cylinder  surface.  When  the 
cylinder  temperature  exceeds  that  of  the  air,  air  adjacent  to  the  cylinder  surface  is  heated. 
This  less  dense  air  rises  allowing  cooler  ambient  air  to  replace  it  (Figure  3). 


Cylinder 


Figure  3.  Vertical  Cylinder  Natural  Convection 


In  all  examples  of  heat  transfer,  temperature  differences  give  rise  to  density 
differences.  Therefore,  we  can  say  that  in  the  presence  of  a  body  force,  the  effects  of 
natural  convection  heat  transfer  are  likely  to  co-exist  with  forced  convection  heat  transfer. 
The  opposite  is  not  true.  The  relative  importance  of  these  two  competing  effects  is 
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governed  by  the  ratio  of  suitably  defined  Grashof  number  to  Reynolds  number  (raised  to 
some  suitable  power).  In  the  limit,  as  this  ratio  tends  to  zero,  the  buoyancy  forces  become 
insignificant.  As  the  ratio  tends  to  infinity,  buoyancy  forces  dominate  the  forced 
convection  effects.  If  the  two  competing  effects  are  comparable  in  magnitude,  mixed  or 
combined  convection  results. 

Analysis  of  mixed  convection  flow  fields  requires  knowledge  of  the  two  limiting 
fluid  flow  cases  (i.e.,  induced  buoyant  forces  and  forced  convection  effects).  From 
Gebhart  et  al.  [Ref  6],  the  direction  and  magnitude  of  the  respective  flow  fields  at  any 
local  position  along  the  surface  of  interest  will  dictate  whether  the  flow  remains  attached, 
separates  or  reverses.  Gebhart  et  al.  present  results  for  differing  geometries,  surface 
heating  and  competing  flow  field  strengths  and  directions. 

2.  Thermoacoustic  Streaming 

Early  works  regarding  the  effect  of  sound  and  vibration  on  heat  transfer,  including 
Lee  and  Richardson  [Ref.  2],  Fand  and  Cheng  [Ref.  7],  and  Fand  and  Kaye  [Ref.  8,  9], 
were  instrumental  in  introducing  a  new  fluid  flow  phenomenon  called  thermoacoustic 
streaming.  All  of  the  above  mentioned  authors’  research  utilized  a  horizontal  heated 
cylinder  in  air  subjected  to  a  horizontal  and  transverse  stationary  sound  field.  The 
boundary  layer  velocity  fields  established  for  such  an  arrangement,  when  not  subjected  to 
a  sound  field,  include  azimuthal  and  radial  velocities  as  expected  in  natural  convection. 
The  former  give  rise  to  natural  convection  currents  tangential  to  the  cylinder  surface 
while  the  latter  satisfy  continuity  requirements  and  are  directed  toward  the  cylinder.  The 
magnitude  of  these  tangential  currents  far  outweigh  the  seemingly  insignificant  radial 
currents.  When  a  sound  field  is  introduced,  a  new  boundary  layer  flow  pattern  is 
superimposed  upon  the  natural  convection  flow.  This  acoustic  streaming  pattern, 
illustrated  by  Fand  [Ref  10],  is  presented  in  Figure  4.  What  appear  to  be  two  distinct 
flow  patterns,  inner  and  outer  streaming  layers,  are  shown.  The  thickness  of  the  inner 
layer  is  exaggerated  for  clarity.  The  order  of  magnitude  of  the  streaming  velocities 
depends  upon  the  intensity  of  the  sound  field  and,  for  much  of  the  above  work. 
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corresponds  to  the  radial  component  of  the  natural  convection  currents.  The  outer 
streaming  flow  aids  the  radial  flow  at  the  top  and  bottom  of  the  cylinder  but  opposes  it  on 
the  sides.  If  the  basic  boundary  layer  flow  pattern  around  this  cylinder  is  changed,  it  is 
reasonable  to  assume  that  the  heat  transfer  may  also  be  affected. 


Figure  4.  Acoustic  Streaming  Pattern 


In  an  isothermal  experiment  this  acoustic  streaming  pattern  behavior  is  well 
documented.  However,  when  applied  to  the  problem  where  temperature  gradients  exist, 
the  streaming  pattern  deviates  from  this  theory.  Several  explanations  are  given  in  Fand 
[Ref.  10].  For  their  specific  experiment,  Fand  and  Kaye  [Ref.  8]  illustrate  pictorially  the 
effects  of  sound  on  the  well  known  convection  flow  pattern.  Two  plumes  or  vortices 
form,  separate  and  grow  above  the  cylinder  as  the  sound  intensity  is  increased.  The 
traditional  natural  convection  boundaiy'  layer  flow  is  disturbed.  Thermoacoustic 
streaming  boundaiy'  layer  flow,  characterized  by  the  two  vortices,  results.  The  streaming 
flows  aid  heat  transfer  at  the  top  and  bottom  but  hinder  it  on  the  sides.  Therefore,  a 
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balance  is  established  whereby  the  overall  heat  transfer  rate  remains  unchanged.  Lee  and 
Richardson’s  [Ref  2]  experimental  results  corroborate  this  theory  noting  that  sound  field 
intensities  below  a  “critical”  sound  pressure  level  (SPL)  have  no  effect  on  the  natural 
convection  heat  transfer  rate.  Once  the  sound  intensity  reaches  the  critical  SPL,  however, 
the  point  of  separation  for  the  thermoacoustic  vortices  has  moved  down  the  sides  of  the 
cylinder.  Streaming  velocities  separate  the  boundary  layer  in  the  horizontal  direction 
after  which  the  buoyant  forces  sweep  the  fluid  upwards  until  the  previous  area  of  local 
reduction  in  heat  transfer  is  eliminated.  The  delicate  heat  transfer  balance  is 
compromised  yielding  an  increase  in  overall  heat  transfer  fi-om  the  surface.  Control  of 
the  heat  transfer  passes  from  pure  natural  convection  to  acoustic  streaming. 

The  above  discussion  and  associated  references  outline  the  theory  and 
experimental  results  regarding  attached  boundary  layer  flow  for  heat  transfer  from  a 
horizontal  cylinder  under  the  influence  of  a  zero-mean  oscillating  flow.  Little  data  exists 
concerning  the  effect  such  a  sound  field  has  on  heat  transfer  from  a  vertical  cylinder. 
More  importantly  for  this  research,  the  theory  and  experimental  data  available  for 
separated  flow  heat  transfer  from  a  vertical  cylinder  exposed  to  a  zero-mean  oscillating 
flow  is  sparse  to  say  the  least. 
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III.  FLUID  MECHANICS 


To  extract  meaningful  results  from  experiments  conducted  by  numerous 
researchers  using  a  variety  of  physical  constants  and  geometric  dimensions,  it  has  become 
common  to  present  them  in  a  general  nondimensional  form.  This  methodology  reduces 
the  number  of  variables  and  links  one  researcher’s  data  to  another,  thereby,  enabling  the 
engineer  the  opportunity  to  efficiently  wdthdraw  relevant  conclusions  regarding  the 
phenomenon  for  design  purposes.  Mapping  the  results  in  terms  of  dimensionless 
parameters  also  gives  the  researcher  insight  into  those  areas  or  boundaries  not  yet 
completely  analyzed.  The  important  parameters  necessaiy  to  explain  the  phenomenon  of 
heat  transfer  from  a  cylinder  in  an  zero-mean  oscillating  flow  are  presented  below. 
Different  fluid  flow  regimes  result.  This  work  attempts  to  extract  heat  transfer 
correlations  from  experimental  data  collected  within  the  incompressible,  separated  flow 
regime. 

A.  RADIAN  WAVELENGTH 

Two  obvious  length  scales  can  be  derived  from  periodic  oscillating  flow.  First, 
the  radian  wavelength  is  defined  as 


In  order  to  maintain  our  assumption  of  incompressibility,  the  radian  wavelength  should  be 
larger  than  the  characteristic  dimension  of  the  wire,  radius,  a.  Wire  radius  is  an  easily 
measured  quantity  as  is  radian  wavelength  when  presented  in  the  form 


A 

In 


2n  0) 


(2) 
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The  dimensionless  parameter  that  results  shall  be  designated,  %.  For  incompressibility  we 
require 


a  acD 


(3) 


B.  DISPLACEMENT  AMPLITUDE 

The  fluid  particles  subjected  to  a  time  varying  sound  field  will  experience  a 
displacement  amplitude,  A .  This  is  the  second  logical  acoustic  length  scale.  One  can 
now  begin  to  appreciate  the  complexity  of  this  problem.  In  addition  to  the  common 
characteristic  geometric  length  scale  (e.g.  wire  radius,  a)  found  in  all  bluff  body  fluid  flow 
experiments,  two  additional  acoustic  length  scales  have  been  introduced. 

In  order  for  flow  near  the  wire  surface  to  remain  attached,  we  expect  minimal 
particle  displacement  relative  to  the  wire  diameter  near  the  wire  surface.  The  amplitude 
parameter,  e,  defined  as 


(4) 


expresses  the  distance  a  fluid  particle  moves  relative  to  the  cylinder  radius.  Its  magnitude 
helps  predict  the  type  of  flow,  separated  or  attached,  along  the  wire  surface.  When  £>  1 
we  anticipate  separated  flow.  For  ease  of  measurement  we  can  expand  this  equation. 
Consider  a  sinusoidal  time  vaiying  signal.  A  direct  relationship  exists  between  particle 
amplitude  and  velocity  as  seen  from  the  following  simple  scenario.  If  particle  amplitude 
is  defined  by 

A  =  Acos{cot)  (5) 

then  particle  velocity  follows  as 
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uJ-^ 

dt 


(6) 


Consequently, 


U  Aasmicot)  (7) 

The  velocity  amplitude,  Uo,  where 


U„  =  Ao) 


(8) 


produces  the  direct  relationship 


CO 


(9) 


One  may  now  measure  velocity  amplitude  along  with  frequency,  or  we  can  carry  the 
scenario  one  step  further.  Morse  and  Ingard  [Ref.  11]  present  a  derivation  for  the 
relationship  between  fluid  velocity  and  pressure  change  for  adiabatic  compression  of  a 
perfect  gas  as 


(10) 


Putting  all  of  this  together  in  a  form  easily  measured  by  experiment,  we  have 


s  = 


c 

aco 


(H) 
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The  amplitude  parameter  by  itself  is  an  indication  of  separated  or  attached  flow. 
Additionally,  when  combined  with  the  Mach  number  results  as 


M-X£  = 


c 


(12) 


I 

We  can  reinforce  the  assumption  of  incompressible  flow  when  M«l.  In  much  of  the 
literature  the  amplitude  parameter  is  also  expressed  as  the  Keulegan-Carpenter  number, 

KC=fC£  (13) 


another  dimensionless  parameter. 

C.  FREQUENCY  PARAMETER 

It  is  important  to  include  the  effects  of  fluid  viscosity  when  discussing  velocity 
and  displacement.  Recall  Stokes’  second  problem.  This  addresses  the  effect  an 
oscillating  plate  has  on  an  adjacent  viscous  fluid.  A  bovmdary  layer  thickness  (or  Stokes 
layer)  of  order 


(14) 


dictates  how  far  into  the  fluid  the  disturbances  created  by  the  oscillating  plate  can  be  felt. 
Solutions  to  Stokes’  second  problem  show  that  at  a  distance  greater  than  approximately 
105  from  the  plate  surface,  the  amplitude  of  fluid  particle  oscillation  has  decayed  to 
within  one  percent  of  the  bulk  fluid  velocity.  For  boundary  layer  flows,  this  Stokes  layer 
thickness  is  of  the  same  order  of  magnitude  as  the  inner  streaming  layer  thickness 
introduced  earlier. 
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A  frequency  parameter  defined  by 


a^o) 

V 


(15) 


relates  the  characteristic  dimension  of  the  body  to  the  Stokes  layer  thickness. 

It  has  been  shown  that  for  A  »1,  a  new,  steady,  acoustic  streaming  velocity 
develops  within  an  outer  layer.  This  steady  flow  component  is  of  order  O(eUo)  and 
appears  as  an  apparent  slip  velocity  on  the  cylinder  surface.  It  decays  to  zero  within  the 
outer  layer.  Reynolds  stresses  linked  to  the  oscillating  viscous  flow  in  the  Stokes  layer 
are  the  driving  force  for  this  new  steady  flow.  Stuart  [Ref  12]  provides  a  mathematical 
proof  of  this.  In  studies  of  oscillatory  flows,  sometimes  the  parameter 


(2a)V 

V 


(16) 


is  used  to  characterize  the  flow. 

D.  REYNOLDS  NUMBER 

Reynolds  number  has  long  been  used  as  a  parameter  to  delineate  flow  regimes 
within  a  given  environment.  Accurate  use  of  this  dimensionless  parameter  requires  the 
proper  definition  of  velocity  and  characteristic  length.  Here,  the  characteristic  length  is 
defined  as  the  cylinder  radius,  a.  And,  in  an  attempt  to  find  appropriate  heat  transfer 
correlations,  two  slightly  different  definitions  for  velocity  have  been  employed.  The  first 
comes  from  the  steady,  acoustic  streaming  velocity  giving  us  the  streaming  Reynolds 
number  as 


K  = 


V 


(17) 
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After  some  manipulation,  the  streaming  Reynolds  number  can  be  related  to  the  acoustic 
parameters  introduced  above  as 


R.  = 


(18) 


The  second  comes  fi-om  the  fluid  particle  oscillation  velocity  resulting  in 


Re 


a 


V 


Re^ 


Uoi2a) 

V 


(19) 


Again,  after  some  rearranging 

Re,  =  ,  Re^  =  =fiKC  (20) 

E.  FLOW  SEPARATION 

Fluid  dynamicists  (Bearman,  Sarpkaya  and  Williamson  to  name  a  few)  have  long 
considered  the  mechanics  of  flow  separation  from  bluff  bodies  in  an  oscillating  cross 
flow.  Active  areas  of  research  today  include:  forces  induced  by  the  separated  flow,  a 
pictorial  representation  of  the  wake  patterns  and  numerical  solutions  that  describe  the 
pointwise  velocities  within  these  wakes  for  both  two  and  three  dimensional  flows. 
Theoretical  hypotheses  and  experimental  results  have  educated  engineers  about  the 
effects  these  flow  patterns  have  on  the  lift  and  drag  forces  imposed  upon  bodies.  The 
practical  application  of  these  studies  has  dealt  mainly  with  the  hydrodynamics  of  marine 
structures  and  vehicles.  Extension  of  these  results  to  heat  transfer  effects  has  been 
minimal. 

In  1980  Honji  [Ref.  13]  presented  a  paper  on  his  observations  of  isothermal  three- 
dimensional  separated  flow  from  a  vertical,  oscillating,  circular  cylinder  in  water.  The 
resulting  photographs  of  “streaked”  flow  from  various  vantage  points  removed  some  of 
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the  mystery  surrounding  the  flow  instabilities  associated  with  large  amplitude  oscillations 
(Figure  5). 


Figure  5.  Isothermal  “Streaked”  Flow 


Williamson  [Ref.  14]  uses  several  vortex  flow  regimes,  each  segregated  by 
Keulegan-Carpenter  number,  to  present  pictorially,  the  formation  and  movement  of 
vortices,  and  graphically,  the  variation  in  lift  and  drag  forces.  He  concludes  that  with  an 
oscillating  flow,  fluid  reversals  have  a  major  effect  on  fluid  induced  forces.  As  oscillation 
amplitude  is  increased,  the  frequency  per  cycle  of  vortex  shedding  increases  and  the  angle 
of  vortex  departure  relative  to  the  incident  wave  varies,  thereby  yielding  higher  induced 
forces  whose  oscillations  are  themselves  periodic.  Furthermore,  the  regularity  of  the 
vortex  pattern  appears  to  break  down  once  a  certain  amplitude  of  oscillation  is  reached, 
thereby  generating  a  discontinuous  and  weakened  wake  pattern.  Williamson 
hypothesizes  that  three  dimensional  effects  along  with  destructive  vonicit}’  between  shed 
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vortices  created  in  the  present  half  cycle  and  those  created  several  cycles  previous  may  be 
the  cause. 

So,  when  does  flow  over  a  bluff  body  separate  and  at  what  Keulegan-Carpenter 
number  can  vortex  shedding  be  observed?  Faltinsen  [Ref  15],  referencing  the  works  of 
Bearman,  Sarpkaya  and  Williamson,  summarizes  the  answers  to  some  of  these  questions. 
It  is  commonly  agreed  upon  that  once  vortices  can  be  seen,  separation  has  occurred.  That 
is  not  to  say,  however,  that  the  vortices  are  necessarily  carried  away  from  the  vicinity  of 
the  body.  At  very  low  Keulegan-Carpenter  numbers  fluid  reversal  generates  mixing  of 
vortices  formed  in  the  previous  half-cycle  with  those  of  the  present.  The  circulation  of 
these  are  of  opposite  sense  and  may  cancel  each  other  out.  Sarpkaya  [Ref  16]  reported 
separation  at  a  Keulegan-Carpenter  number  of  1.25.  Williamson  [Ref  14]  has 
demonstrated  vortex  shedding  at  Keulegan-Carpenter  numbers  as  low  as  four. 

It  is  clear  that  significant  work  has  been  documented  regarding  the  fluid 
mechanics  associated  with  flow  separation  and  vortex  shedding  over  a  bluff  body.  It  is 
logically  hypothesized  here,  that  similar  flow  instabilities  will  arise  when  temperature 
gradients,  in  addition  to  large  amplitude  fluid  oscillations,  are  imposed  upon  a  fluid 
medium  near  a  heated  cylinder  surface.  For  this  study  then,  it  would  seem  safe  to  say 
that  in  addition  to  the  previously  mentioned  requirements  for  separation,  if  Keulegan- 
Carpenter  numbers  above  four  are  observed,  vortex  shedding  occurs. 


IV.  EXPERIMENT 


A.  PURPOSE 

Not  all  problems  of  science  can  be  solved  analjlically.  Many  times  it  is  through 
experiment  that  the  true  nature  of  a  scientific  phenomenon  is  best  explained  and 
xmderstood.  A  small  sample  of  the  considerable  theory  and  experimental  data  concerning 
the  effects  of  sound  on  attached  boundary  layer  flow  heat  transfer  from  a  horizontal 
cylinder  has  been  presented.  The  research  in  this  area  is  hardly  exhaustive,  however,  it 
far  outweighs  the  available  material  detailing  separated  flow  heat  transfer  from  a  vertical 
cylinder  in  the  presence  of  a  zero-mean  oscillating  sound  field.  What  is  available 
concerning  separated  flows  in  an  isothermal  environment  has  also  been  presented  briefly 
to  introduce  the  relationship  between  oscillating  fluid  amplitude  and  vortex  shedding. 

Similar  geometric,  acoustic  and  fluid  mechanics  parameters,  presented  above,  are 
useful  in  quantifying  experimental  results  in  both  cases.  Such  a  standardization  allows 
effective  comparisons  between  different  experimental  works  and  provides  researchers  a 
“map”  on  which  to  track  the  evolution  of  analyses  xmtil  a  complete  imderstanding  is  at 
hand.  It  is  at  this  point  that  a  confident  assessment  of  the  usefulness  of  thermoacoustics 
in  commercial  heat  exchanger  design  can  be  made. 

Analysis  of  mean  external  flows  over  bluff  bodies  leads  one  to  solve  the  boundary 
layer  equations.  Friction  and  convective  heat  transfer  correlations  developed  from  these 
solutions  have  proven  invaluable  to  modem  engineers  in  their  design  of  new  technologies 
ranging  from  spacecraft  to  golf  balls.  On  the  other  hand,  little  has  been  attempted 
theoretically  or  experimentally,  to  analyze  the  heat  transfer  behavior  from  bodies  in  zero- 
mean  oscillating  flow  fields.  Simple  analogies  from  mean  flow  behavior  are  no  longer 
directly  applicable. 

Fand  and  Kaye  [Ref  9]  describe  related  work  done  by  Kubanskii  in  1952.  His 
experimental  arrangement  placed  the  stationary  sound  field  horizontal  and  parallel,  vice 
transverse,  to  the  axis  of  a  horizontal  heated  cylinder.  He  deduced  that  fluid  near  the 
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cylinder  surface  moved  from  the  antinodes  to  the  nodes  at  which  point  the  flows  meeting 
from  opposite  directions  collided  to  form  a  single  stream  that  moved  away  perpendicular 
to  the  cylinder  in  a  vortex-type  or  separated  flow  fashion.  These  experiments  aim  to 
create  similar  vortex-type  flow  but  from  a  vertical  cylinder  subjected  to  a  transverse 
sound  field.  The  anticipated  driving  force  for  these  experiments  is  forced  convection 
rather  than  the  free  convection  found  by  Kubanskii. 

B.  EQUIPMENT 

The  approach  used  here  will  be  similar  to  that  used  by  Harder  [Ref  17].  Much  of 
the  supporting  equipment  is  the  same,  with  the  biggest  difference  coming  in  the  diameter 
of  the  test  cylinder.  A  slightly  modified  acoustic  chamber  and  electric  circuit  enable  the 
smooth  integration  of  the  new  test  cylinder  into  the  experimental  setup. 

1.  Test  Cylinder 

Recent  heat  transfer  correlations  have  been  developed  using  small  diameter 
vertical  cylinders  in  a  stationary  sound  field.  George  Mozerkewich  [Ref  18]  conducted 
transient  heat  transfer  rate  analysis  on  rhenium  and  titanium  wires  of  diameter  ranging 
from  0.125  to  2  mm.  Don  Harder  [Ref  17]  utilized  a  5  mm  diameter  copper  cylinder 
with  an  imbedded  cartridge  heater  to  develop  steady  state  heat  transfer  rate  solutions.  In 
the  present  research,  a  0.05  mm  (and  0.127  mm)  diameter  platinum  wire  acts  as  both  the 
heater  and  thermometer  for  steady  state  heat  transfer  analysis.  One  of  the  more  attractive 
aspects  of  platinum  is  its  nearly  linear  temperature-resistance  relation.  Upon  supplying  a 
known  current  through  the  approximately  75  mm  length  of  wire,  the  voltage  drop  can  be 
measured  and  from  Ohm’s  Law  resistance  is  easily  calculated. 

The  wire  was  attached  at  each  end  to  a  copper  stud.  A  small  indentation  was 
made  on  the  flat  face  of  each  copper  stud  to  allow  pooling  of  the  flux  and  platinum  for 
increased  mechanical  integrity.  The  stud  was  then  washed  in  a  nitric  acid  solution.  This 
provided  a  “sticky”  surface  onto  which  the  flux  could  better  adhere.  The  platinum  wire 
was  then  soldered  at  either  end  to  the  threaded  copper  studs.  In  order  to  maintain  a 
relatively  constant  tension  in  the  wire  during  both  calibration  and  testing,  the  lower 
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copper  stud  remained  free  to  move  vertically,  thereby,  permitting  very  short  extension 
and  contraction  movement  during  temperature  changes. 

Imposing  a  temperature  versus  resistance  linear  regression  model  upon  the 
calibration  data  (App.  A),  vvire  temperature  is  determined  for  all  experimental  runs  from 
the  calculated  resistance.  The  importance  of  accurate  resistance  measurements  across  the 
platinum  wire  cannot  be  overemphasized.  Very  low  resistance  copper  and  brass  circuit 
connections  were  used  both  in  calibration  and  testing.  A  Kikusui  Model  PAR  160A 
regulated  DC  power  supply  provided  circuit  voltage  and  current  control.  Voltage  and 
current  resolution  are  approximately  5  mV  and  5  mA,  respectively. 

2.  Acoustic  Chamber 

Current  development  of  thermoacoustic  refrigerators  primarily  use  either  U- 
shaped  or  straight  circular  cross-section  housings  to  support  the  acoustic  driver,  heat 
exchangers  and  stack.  Sound  is  confined  quite  well,  and  standing  waves  can  be  easily 
created  in  such  simple  structures.  Figure  6  depicts  a  sketch  of  the  acoustic  chamber  used 
in  these  experiments.  Harder  [Ref.  17]  aptly  describes  the  chamber’s  physical 
characteristics  prior  to  a  few  design  changes  all  of  which  are  described  below. 

Using  approximate  dimensions,  a  1.8  m  plexiglass  tube  of  9  cm  outside  diameter 
and  6  mm  wall  thickness  is  capped  at  one  end  by  a  JBL  Model  2490H  acoustic  driver  and 
a  moveable  end  plate  at  the  other.  The  tube  is  supported  by  stanchions  at  various 
positions  along  its  length.  An  Endevco  piezoresistive  pressure  transducer  embedded  in 
the  moveable  endplate  provides  data  for  pressure  ratio  and  sound  pressure  level 
determination  within  the  chamber.  An  o-ring  seal  sunounding  this  endplate  permits  its 
easy  movement  into  and  out  of  the  chamber  while  confining  the  sound  wave.  Along  the 
length  of  the  chamber  are  bored  two  sets  of  diagonally  opposed  12  mm  diameter  holes. 
During  testing,  the  platinum  wire  assembly  is  placed  at  one  of  these  locations  while  the 
other  is  blank-capped  except  for  ambient  temperature  measurements  taken  inside  the 
chamber. 
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Figure  6.  Acoustic  Chamber  and  Accessories 


3.  Electronics 

Several  pieces  of  electronics  equipment  are  used  to  generate  the  acoustic  signal 
and  then  to  monitor  some  of  its  properties.  Figure  7  depicts  the  intercoimecting 
relationships  amongst  this  equipment.  Beginning  with  a  Hewlett  Packard  Model  33 120 A 
waveform  generator,  the  desired  waveform  (sinusoidal),  frequency  and  amplitude  are 
selected.  This  signal  is  fed  to  a  Techron  Model  7540  power  supply  amplifier  to  boost  the 
signal  amplitude  imtil  a  desired  pressure  ratio  or  sound  pressure  level  (SPL)  is  reached. 
The  amplified  waveform  enters  the  acoustic  chamber  via  a  JBL  Model  2490H  acoustic 
driver  located  at  one  end  of  the  chamber.  A  pressure  antinode  is  established  at  the 
moveable  endplate  located  at  the  opposite  end  of  the  acoustic  chamber.  The  standing 
wave  pressure  signal  is  sensed  by  the  pressure  transducer  and  amplified  100  times  before 
it  is  displayed  on  a  Hewlett  Packard  Model  34401 A  multimeter  as  an  AC  voltage.  It  is  at 
this  point  that  the  pressure  ratio  and  SPL  within  the  chamber  can  be  determined. 
Simultaneously,  the  amplified  pressure  signal  is  viewed  on  an  oscilloscope  to  monitor 
gross  waveform  regularity.  In  series  with  the  oscilloscope,  a  Hewlett  Packard  Model 
3562A  dynamic  signal  analyzer  can  record  the  signal  frequency  and  amplitude  for  all 
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harmonics.  An  indication  of  harmonic  interference  from  the  signal  analyzer  helps 
determine  exactly  which  frequency-amplitude  combinations  establish  the  optimum 
acoustic  conditions  in  the  chamber  prior  to  the  steady  state  heat  transfer  rate  experimental 
run. 


Figure  7.  Acoustic  Electronics  Package 


Once  the  desired  acoustic  conditions  are  established,  the  following  additional 
equipment  (Figure  8)  is  used  to  establish  and  monitor  the  heat  transfer  characteristics  of 
the  experiment.  Ambient  temperature  sensing  within  the  acoustic  chamber  is  performed 
by  manually  inserting  a  T-type  thermocouple  into  the  center  of  the  chamber  through  an 
upstream,  normally  capped,  hole  at  the  appropriate  time.  A  Keithley  Model  740  system 
scanning  thermometer  then  measures  the  temperature  to  an  accuracy  of  ±0.5°C. 

The  electrical  circuit  designed  to  provide  the  desired  temperature  versus  resistance 
characteristics  for  the  platinum  wire  consists  of  a  known  fixed  resistance  and  the 
platinum  wire  connected  in  series  with  a  power  supply.  The  resistor  is  rated  at  98.3  mQ 
with  an  accuracy  of  0.026  %  and  acts  both  to  minimize  the  power  drawn  by  the  platinum 
wire  and  to  provide  a  means  for  accurate  circuit  current  measurement.  A  Hewlett 
Packard  Model  3478A  digital  multimeter  connected  in  parallel  across  both  the  resistor 
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Figure  8.  Heat  Transfer  Electronics  Package 
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and  wire  accurately  measures  voltages  up  to  a  5  '/2  digit  resolution.  These  measurements 
are  subsequently  used  to  calculate  circuit  current  and  wire  resistance. 

C.  PROCEDURE 

Prior  to  inserting  the  platinum  wire  into  the  acoustic  chamber,  a  series  of  tests 
were  conducted  to  identify  the  resonant  frequencies  for  various  chamber  end  lengths,  L. 
With  two  sets  of  diagonally  opposed  holes  bored  into  the  chamber  wall,  the  second  set 
being  50  cm  offset  from  the  first  (Figure  6),  end  lengths  were  incremented  by  1  cm  from  a 
maximum  of  73  cm  (123  cm)  to  a  minimum  of  42  cm  (92  cm)  by  adjusting  the  position  of 
the  moveable  endplate.  Sine  waves  starting  at  300  Hz,  the  lower  end  of  the  acoustic 
driver  operating  range,  were  introduced  into  the  chamber.  Frequency  was  raised  until  a 
resonant  condition  existed  as  evidenced  by  the  wave  shape  seen  on  the  oscilloscope  and 
peak  pressure  transducer  output  as  read  from  a  multimeter  in  root  mean  square  (RMS) 
volts  AC.  This  resonant  conditions  occurred  at  roughly  100  Hz  increments  up  to 
approximately  1900  Hz,  the  upper  end  of  the  acoustic  driver  operating  range.  For  each 
resonant  state,  the  frequency,  end  length,  peak  RMS  voltage  for  the  fundamental  and 
second  harmonic  frequencies,  and  maximum  pressure  transducer  output  were  recorded. 
The  strength  of  the  second  harmonic  provided  a  measure  of  harmonic  interference  at  the 
velocity  antinode.  This  information  established  which  frequency  and  end  length 
combinations  to  use  for  actual  testing. 

The  geometry  of  the  experimental  setup  plays  an  important  role  in  proper  test 
wire  positioning.  For  a  given  resonant  frequency  and  end  length  combination,  the 
standing  wave  velocity  created  can  be  depicted  as  in  Figure  9.  Test  wire  placement  at 
maximum  sound  velocity  (antinode)  should  produce  the  greatest  heat  treinsfer  from  the 
wire.  Here  a  simple  analogy  to  forced  convection  may  be  helpful.  If  instead  of  an 
acoustic  driver,  a  fan  were  placed  at  one  end  of  the  chamber  (open-ended  at  the  opposite 
end  in  this  case),  the  speed  of  the  fan  would  effect  the  heat  transfer  rate.  The  higher  the 
fan  speed,  and  therefore  air  velocity  in  the  vicinity  of  the  wire,  the  higher  would  be  the 
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heat  transfer  rate.  Similarly  for  thermoacoustics,  the  greater  the  medium  velocity  near  the 
wire,  the  greater  the  anticipated  heat  transfer  rate. 

As  is  evident  from  Figure  9  any  odd  multiple  of  quarter  wavelength  (i.e.,  A,/4, 
3X-/4,  5X,/4  etc.)  lies  at  a  velocity  antinode.  From  the  equation 

1  =  —;  /j  =  U5...  (21) 

4 


we  can  determine  whether  or  not  a  given  resonant  condition  places  the  test  wire  at  a 
velocity  antinode.  Rearranging,  we  get 


«  =  1,3,5... 


(22) 


Now,  knowing  which  combinations  of  resonant  frequency  and  end  length  were  best  suited 
for  the  fixed  geometric  conditions,  the  maximum  signal  amplitude  yielding  minimal 
harmonic  interference  can  be  found. 


Figure  9.  Standing  Wave  Velocity  Profile  in  a  Circular  Cylinder 


P.  D.  Richardson  [Ref  19]  points  out  that  rather  than  measuring  velocity 
amplitude  with  a  hot-wire,  it  is  more  accurate  and  convenient  to  measure  pressure 
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amplitude  with  a  microphone,  as  done  here,  since  the  velocity  is  directly  proportional  to 
the  pressure  ratio.  Recalling  that  the  resonant  sound  wave  pressure  is  90  degrees  out  of 
phase  from  the  velocity,  a  pressure  antinode  exists  at  the  moveable  endplate  where  it  is 
detected  by  the  microphone.  The  microphone  output  is  amplified  100  times.  The  signal 
is  then  measured  by  a  Hewlett  Packard  34401 A  multimeter  as  RMS  volts  AC.  The 
following  equation 


gives  the  maximum  pressure  within  the  chamber.  The  amplified  pressure  transducer 
sensitivity,  S,  is  equal  to  0.738  mV/Pa.  Solving  for  pressure,  we  get 


Using  101 .35  kPa  for  mean  ambient  pressure.  Pm,  we  can  find  the  pressure  ratio  from 


PR  =  ^*m% 


(25) 


and  sound  pressure  level  (SPL)  from 


SPL  =  201og,o 


(26) 
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where  Pref  is  the  reference  pressure,  taken  to  be  20  pPa  for  gases.  In  order  to  cover  a 
wide  range  of  fluid  particle  amplitudes,  the  pressure  ratios  used  during  the  experiment 
varied  from  0.9  to  2.9  percent  in  0.1  percent  increments  or  up  to  the  point  where 
harmonic  interference  became  greater  than  approximately  ten  percent. 

Once  the  optimum  acoustic  setpoints  were  found,  the  calibrated  test  wire  was 
positioned  within  the  acoustic  chamber.  Movement  of  the  wire  from  its  calibration 
chamber  to  the  experimental  acoustic  chamber  required  very  delicate  maneuvering.  Any 
excessive  torsional  or  axial  stress  exerted  on  the  copper  to  platinum  wire  solder  joint 
separated  the  joint.  After  the  vvire  was  positioned,  the  complete  electrical  circuit  was 
established  including  all  measuring  instrumentation.  At  this  point,  a  comprehensive  plan 
was  set  forth  from  which  repeated  experimental  data  runs  could  be  made  with  minimal 
procedural  variance  between  them.  The  following  description  details  one  such 
experimental  data  run. 

After  an  initial  equipment  warm  up  period,  one  frequency  and  end  length 
combination  was  set.  Acoustic  signal  amplification  was  adjusted  until  the  desired 
pressure  ratio  was  obtained.  Fine  tuning  the  frequency  ensured  a  resonant  condition 
within  the  acoustic  chamber.  The  signal  characteristics  were  monitored  on  both  the 
oscilloscope  and  frequency  analyzer  in  order  to  check  the  level  of  harmonic  interference. 
Signal  amplifier  settings  were  noted  and  then  the  signal  was  secured.  After  logging  an 
ambient  temperature  reading  within  the  chamber,  a  preset  source  current  and  voltage  were 
introduced  into  the  circuit.  As  the  wire  heated  up,  the  previously  marked  signal  amplifier 
settings  were  re-established.  Once  steady,  after  apjDroximately  15  seconds,  voltage 
measurements  across  the  series  precision  resistor  and  test  wire  were  taken.  As  soon  as 
these  readings  were  recorded,  the  power  circuit  source  and  signal  amplifier  were  secured. 
A  second  ambient  temperature  reading  within  the  chamber  was  logged.  Differences 
between  initial  and  final  ambient  temperature  readings  could  be  as  high  as  0.3  degrees. 
All  of  the  above  readings  were  then  entered  into  a  spreadsheet  where  acoustic  and  heat 
transfer  parameters  were  calculated.  Through  trial  and  error  the  above  procedure  was 
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repeated  until  the  correct  source  current  and  voltage  were  found  for  the  desired 
temperature  difference  between  the  ambient  and  wire  surface.  Temperature  differences  of 
8,  16,  and  24  degrees  (±0.5  degrees)  were  targets  for  each  pressure  ratio.  At  some  of  the 
highest  pressure  ratios  for  the  127  pm  diameter  wire,  temperature  differences  of  8,  12, 
and  16  degrees  were  used  because  measured  wire  voltage  fluctuations  were  unacceptably 
high  at  24  degrees.  In  addition  to  three  different  temperature  runs  per  pressure  ratio, 
each  voltage  reading  within  a  run  was  taken  in  triplicate.  Averaging  the  results  from  three 
measurements  helped  minimize  the  error  for  a  given  parameter  within  a  data  run. 

The  above  procedure  was  used  at  every  pressure  ratio  and  for  every  frequency  and 
end  length  combination  selected.  In  addition,  nearly  the  same  procedure,  minus  the  sound 
signal,  was  followed  to  determine  the  heat  transfer  coefficient  for  natural  convection. 
The  results  were  compared  to  those  presented  by  Gebhart  et  al.  [Ref  6]  for  natural 
convection  from  a  semi-infinite  thin  vertical  wire  in  air  to  validate  the  procedure.  Our 
enclosed  container  geometry  does  not  accurately  replicate  previous  experimental  setups, 
therefore,  our  results  do  not  match  those  presented  by  Gebhart.  However,  our  results 
were  of  the  same  order  of  magnitude,  This  offered  some  validity  to  our  experimental 
setup  and  procedure. 
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V.  RESULTS  AND  DISCUSSION 


This  study  investigates  the  heat  transfer  characteristics  of  a  thin  vertical  cylinder 
in  an  air  medium  subjected  to  a  high  amplitude  zero-mean  oscillating  flow.  Two 
different  wire  diameters  were  individually  tested,  50  pm  and  127  pm.  Both  wires  were 
exposed  to  pressure  ratios  incremented  from  0.9  to  2.9  %  while  ensuring  minimal 
harmonic  interference.  Additionally,  at  each  pressure  ratio  setting,  data  for  three  separate 
temperature  differences  between  the  ambient  and  wire  surface  were  taken.  In  all,  nearly 
400  data  points  were  collected.  Experimental  measurements,  applicable  heat  transfer  and 
acoustic  parameter  calculations  and  certain  parameter  uncertainties  are  presented  in  a 
spreadsheet  format  foimd  in  Appendix  D. 

Use  of  such  thin  wires  all  but  guaranteed  separated  flows  in  the  high  amplitude 
oscillation  environment.  The  calculated  values  for  the  amplitude  parameter,  (i.e., 
Keulegan-Carpenter  number),  which  were  of  0(10),  best  corroborated  this  supposition. 
An  early  attempt  was  made  to  use  a  relatively  thick  254  pm  diameter  wire  with  little 
success.  With  this  wire,  resistance  measurements  were  on  the  order  of  100  mQ  with  a 
1°C  temperature  change.  This  corresponded  to  a  0.6  mQ  resistance  change.  Such  small 
values  for  resistance  generated  unacceptably  high  uncertainties  for  our  measurement 
technique  and  available  equipment,  thereby  forcing  the  use  of  thinner  wires. 

The  first  useful  experiments  were  conducted  on  a  50  pm  diameter  wire.  Initially, 
a  frequency  of  1213  Hz  with  a  corresponding  end  length  of  64  cm  was  used.  This 
combination  most  closely  placed  the  maximum  velocity  at  the  fixed  test  wire  position. 
During  testing,  this  combination  also  showed  the  least  harmonic  distortion  up  to  a 
maximum  pressure  ratio  of  2.9  %.  No  other  combination  withstood  the  distortion 
requirements  at  such  a  high  pressure  ratio. 

Next,  the  effect  of  acoustic  chamber  end  length  was  examined.  A  frequency  near 
1213  Hz  but  with  an  end  length  different  than  64  cm  was  desired.  The  one  combination 
available  that  offered  minimal  harmonic  distortion  was  a  frequency  of  1213  Hz  at  an  end 
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length  of  50  cm.  The  calculated  acoustic  and  heat  toinsfer  parameters  from  these  two 
data  sets  were  compared  and  noted  to  be  nearly  identical  for  like  pressure  ratios  and 
temperature  differences.  It  was  concluded  that  different  end  lengths  for  a  given 
frequency,  while  maintaining  a  zero-mean  oscillating  flow  within  the  acoustic  chamber, 
do  not  alter  the  resulting  heat  transfer  correlations.  In  order  to  gather  a  more  diverse  set 
of  data,  two  additional  frequencies  and  associated  end  lengths  were  tested.  The  results  at 
1053  Hz  (L  =  58  cm)  and  564  Hz  (L  =  47  cm)  are  presented  in  Appendix  D. 

Once  the  50  pm  diameter  wire  data  was  organized,  a  systematic  interpretation 
began.  To  examine  the  effects  of  acoustics  on  heat  transfer,  several  Nusselt  number 
correlations  were  investigated  using  the  amplitude  and  frequency  parameters  individually 
and  in  combination  with  each  other.  Nusselt  number  as  a  function  of  Reynolds  number, 
not  streaming  Reynolds  number,  provided  the  best  fit  to  the  raw  data  (Figure  10).  From 
this  data,  a  curve  fit  yielded  the  correlation 

=  0.062  Re'"  (27) 

Nusselt  number  errors  reached  as  high  as  13  %,  while  the  Reynolds  number  error  was 
approximately  5  %.  Again,  the  forced  convection  analogy  can  be  invoked  here  to  help 
explain  these  results.  Reynolds  number  is  a  function  of  particle  velocity  for  a  given 
cylinder  diameter.  It  has  been  shown  that  particle  velocity  is  directly  proportional  to  the 
pressure  ratio,  therefore,  as  the  pressure  ratio  increases  so  does  the  Reynolds  number.  For 
increasing  fluid  velocity  we  expect  greater  heat  transfer  and  that  is  precisely  what  we  see 
in  these  results. 

To  extend  the  range  of  data  and  to  ensure  a  suitable  Nusselt  number  correlation 
had  been  found,  an  additional  set  of  data,  using  the  same  physical  and  acoustic  parameters 
as  above,  was  generated  for  a  127  pm  diameter  wire.  A  graph  of  this  data  for  Nusselt 
number  versus  Reynolds  number  is  shown  in  Figure  1 1 .  For  this  wire,  Nusselt  number 
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errors  were  in  the  range  of  4-15  %,  while  the  Reynolds  number  error  was  approximately  2 

%. 

A  graph  of  the  complete  set  of  data  from  both  wires  is  presented  in  Figure  12.  It 
was  expected  that  the  thicker  wire  data  would  fall  on  top  of  that  for  the  thinner  wire  over 
the  common  Reynolds  numbers,  and  that  the  trend  of  monotomically  increasing  Nusselt 
number  with  Reynolds  number  would  continue  into  a  new  higher  Reynolds  number 
regime.  Neither  of  these  hypotheses  were  accurately  realized.  However,  within  the 
Reynolds  nxunber  overlap  region  for  the  two  wires  one  can  see  that  ev^n  though  the  data 
does  not  strictly  overlap,  the  slopes  of  the  two  data  sets  are  comparable. 

Obviously,  the  thicker  wire  data  yielded  much  more  interesting  results  (Figure 
1 1).  A  steady  increase  followed  by  a  peak  and  subsequent  drop  in  Nusselt  number  occurs 
for  all  three  frequencies.  As  the  frequency  increases,  it  appears  that  the  rate  at  which  the 
Nusselt  number  drops  and  the  magnitude  of  the  drop  also  increase.  Similarly,  for 
increasing  frequency,  the  drop  in  Nusselt  number  occurs  at  a  lower  Reynolds  munber.  A 
discussion  of  these  results  follows. 

For  a  given  wire  diameter  ( i.e.,  127  pm),  increasing  frequency  corresponds  to  an 
increase  in  p.  Physically  speaking,  the  ratio  of  wire  radius,  a,  to  Stokes  layer  thickness,  5, 
increases.  Similarly  stated,  the  effects  of  the  Stokes  fluid  oscillations  do  not  extend  as  far 
into  the  bulk  fluid.  Consequently,  one  could  conclude  that  the  effects  of  Stokes  flow 
interaction  with  the  shed  vortices  is  reduced.  Any  mixing  effect  associated  with  this 
interaction  fades  resulting  in  a  lower  heat  transfer  rate.  In  Figure  11,  three  curves  are 
evident,  one  for  each  of  three  frequencies  with  P  of  0(1).  The  high  frequency  (P)  curve 
sits  below  the  low  frequency  (P)  curve  indicating  less  heat  transfer  for  a  higher  frequency 
as  theorized.  Furthermore,  rather  than  applying  simply  a  Reynolds  number  dependence 
for  the  heat  transfer  correlation  (Equation  27),  we  now  see  a  possible  p  dependence  as 
well.  Reviewing  Figure  10,  the  50.8  pm  diameter  wire  p  dependence  is  not  as  evident. 
Here  p  is  of  0(0.1).  Quite  possibly,  P  becomes  influential  only  when  of  0(1)  or  higher. 
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To  explain  the  sudden  Nusselt  number  drops  at  high  Reynolds  number, 
Williamson’s  results  [Ref.  14]  may  be  helpful.  Recalling  the  work  of  Williamson,  flows 
yielding  high  Keulegan-Carpenter  numbers  evoke  separation  and  varying  vortex  shedding 
wake  patterns  which  he  groups  into  KC  number  regimes.  This  study  includes  KC 
numbers  ranging  from  14  to  200.  If  these  vortex  patterns  proceed  with  some  regularity  a 
steady  increase  in  Nusselt  number  with  KC  may  be  expected.  However,  if  between 
acoustic  half-cycles,  vortex  interference,  out  of  phase  vortex  shedding  along  the  cylinder 
length  or  a  loss  of  regularity  in  vortex  shedding  occur,  one  might  expect  a  disruption  in 
the  steady  heat  transfer  from  the  cylinder.  One  possible  cause,  as  noted  by  Williamson, 
could  be  an  intermittent  and  random  change  in  the  vortex  shedding  mode  (i.e.,  one  side  of 
the  cylinder  suddenly  becoming  the  preferred  side  for  vortex  shedding).  Once  this  change 
occurs  a  regularity  to  the  vortex  shedding  should  reappear.  Figure  1 1  gives  a  limited 
glimpse  at  this  possibility.  The  1213  Hz  data  shows  that  a  minimum  Nusselt  number  is 
reached  after  which  an  upturn  in  the  curve  may  be  indicative  of  a  return  to  vortex 
shedding  regularity.  The  data  to  support  this  is  far  from  conclusive,  but  the  hypothesis  is 
in  keeping  with  Williamson’s  results  regarding  lift  force  coefficient  fluctuations. 
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NUj  vs.  Re^ 
(d  =  127  ^im) 


Figure  11.  Nusselt  Number  versus  Reynolds  Number  (d  =  127  iim) 
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Figure  12.  Nusselt  Number  versus  Reynolds  Number  (both  wires) 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  motivation  for  this  study  was  two-fold.  First  was  the  desire  to  gain  an 
understanding  of  the  underlying  fluid  mechanics  for  oscillating  flow  over  a  bluff  body. 
Second  was  the  need  to  develop  useful  heat  transfer  correlations  for  high  amplitude 
oscillating  flow  over  a  vertical  cylinder  for  application  in  the  design  of  more  efficient 
thermoacoustic  engines.  The  experiments  performed  utilized  a  varying  set  of  acoustic 
parameters  for  two  different  cylinder  diameters.  Two  regimes  of  varying  heat  transfer 
mechanisms  are  discussed. 

Conventional  studies  of  oscillating  flow  over  bluff  bodies  deal  with  large  P  values 
for  varying  Keulegan-Carpenter  numbers.  This  study  has  experimented  with  very  low  p 
and,  subsequently,  two  regimes  of  differing  heat  transfer  mechanisms  have  been 
identified. 

For  p  of  0(0.1)  and  over  a  wide  range  of  moderate  to  high  KC  numbers  (35  to 
200),  Equation  (27)  provides  a  good  curve  fit  to  the  data  and  a  correlation  for  the  Nusselt 
number  as  a  function  of  Reynolds  number  (Figure  13).  In  this  case,  the  oscillating  Stokes 
flow  interacts  with  vortex  shedding  for  a  combined  mechanism  of  heat  transfer.  What 
has  not  been  considered  in  these  experiments  is  the  effect  of  Prandtl  number.  Generally 
speaking,  Nusselt  number  is  a  function  of  both  Reynolds  number  and  Prandtl  number  for 
forced  flow  over  a  cylinder.  This  is  an  area  for  expanded  research. 

For  p  of  0(1)  and  over  a  wide  range  of  moderate  to  high  KC  numbers  (14  to  86), 
no  definitive  correlation  has  been  developed.  An  introduction  to  vortex  shedding  and  its 
effect  on  heat  transfer  in  this  regime  has  been  discussed,  but  additional  data  would  help 
solidify  these  results.  In  particular,  doubling  the  diameter  of  the  test  cylinder,  while 
maintaining  KC  high  enough  for  vortex  shedding,  would  generate  a  set  of  data  by  which 
several  hypothesis  presented  above  could  be  validated.  Would  a  thicker  wire,  and 
therefore  a  larger  p,  generate  curves  with  similar  slope  but  lower  heat  transfer?  If  a  peak 
and  sudden  drop  in  Nusselt  number  occurred,  would  it  be  at  a  lower  Reynolds  number  as 
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the  trend  from  these  results  indicates?  The  answers  to  these  questions  and  validation  of 
the  results  for  p  of  0(1)  will  come  from  future  research. 

Once  a  confident  explanation  of  the  mechanisms  that  control  heat  transfer  for  low 
P  values  has  been  found,  expanded  research  should  include  arrays  of  cylinders.  Analysis 
of  two  or  more  cylinders  of  variable  separation  distance  arranged  as  either  aligned  with 
the  soimd  wave,  staggered  or  side-by-side  is  a  logical  next  step  to  developing  heat 
transfer  correlations  for  operational  heat  exchangers.  A  removable  section  of  the  acoustic 
chamber  would  facilitate  test  platform  interchangeability,  make  calibration  easier  and 
minimize  handling  of  these  delicate  structures.  And  finally,  use  of  an  automated  data 
acquisition  system  would  provide  more  uniformity  in  data  taking  and  speed  up  the  entire 
testing  process. 
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•  raw  data 
—  Linear  Regression 


T - 1 - 1 1 — I — r-|—| - “1 - 1 - 1 - 1 r 

10 

Re^j 


Figure  13.  Curve  Fit  for  Nusselt  Number  versus  Reynolds  Number  (d  =  50  |j.m) 
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APPENDIX  A:  CALIBRATION 


One  of  the  most  basic  measurements  required  in  order  to  effectively  develop  heat 
transfer  correlations  is  temperature.  The  two  temperature  measurements  of  concern  here 
are  of  the  ambient  air  within  the  acoustic  chamber  and  the  surface  of  the  platinum  wire. 
Thermoresistive  elements,  like  a  platinum  wire,  provide  a  nearly  linear  temperature- 
resistance  relation.  A  preliminary  step  to  the  experimentation  was  to  determine  this 
relation  through  calibration.  Ambient  air  temperature  measurements  were  made  with  a 
calibrated  T-type  thermocouple. 

The  platinum  wire  was  mounted  in  a  calibration  chamber  cut  from  the  same 
plexiglass  tube  stock  used  to  make  the  acoustic  chamber.  All  connectors,  electric 
circuitry  and  equipment  were  identical  to  those  used  throughout  experimentation.  Once 
assembled,  as  in  Figure  14,  the  calibration  chamber,  along  with  a  T-type  thermocouple, 
was  immersed  in  a  Rosemont  Engineering  Model  91 3 A  ethyl  glycol  bath.  An  associated 
Rosemont  Engineering  Model  920A  commutating  bridge  provided  the  reference 
temperature  as  measured  by  a  precision  temperature  probe.  A  four  wire  bridge  electrical 
circuit  was  used  to  keep  the  lead  resistance  as  low  as  possible  relative  to  the  platinum 
wire  resistance  [Ref.  20]. 

Numerous  resistance  versus  temperature  measurements  were  taken  between  20°C 
and  50°C,  the  experimental  temperature  range  of  interest.  Tables  1  and  2  show  the 
measured  wire  resistance  and  associated  precision  temperature  values  for  the  two  wire 
diameters.  Column  T1  gives  the  calculated  temperature  resulting  from  the  linear 
regression  equations.  The  error  is  simply  the  percent  difference  between  the  measured 
and  calculated  temperatures.  Figures  15  and  16  show  the  statistical  linear  regression 
curves  superimposed  upon  the  raw  data.  The  following  two  equations  were  used 
throughout  the  data  reduction  phase  of  the  experiment  for  the  50.8  pm  and  127  pm 
diameter  wires  cases,  respectively 

Temperature  (°C)  =  59.8 1 56*Resistance  (Q)  -  224.2799  (A.  1 ) 
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Temperature  (°C)  =  436.1 12’'‘Resistance  (O)  -  261.8139  (A.2) 

Finally,  a  calibration  data  run  was  made  for  the  127  pm  diameter  wire  using  the 
actual  electric  circuit  employed  in  the  experiment.  Assuming  a  very  high  forced 
convection  heat  transfer  coefficient  within  the  ethyl  glycol  bath,  and  applying  a  low 
power  to  the  platinum  wire,  temperature  differences  between  the  wire  surface  and  the 
fluid  are  assumed  insignificant.  Voltage  measurements  were  made  and  wire  resistances 
were  calculated.  Wire  resistance  was  plotted  against  the  precision  temperature  to  generate 
an  independent  calibration  curve  as  seen  in  Figure  16.  The  results  from  these  two 
different  calibration  techniques  were  compared  and  showed  a  small  change  in  slope  and 
intercept  of  the  calibration  curve.  The  four  wire  bridge  calibration  results  from  both 
platinum  wires  were  used  in  all  experimental  calculations  in  order  to  maintain  some 
measure  of  consistency. 


Finally,  calibration  of  the  T-type  thermocouple,  using  the  same  Rosemont 
Engineering  calibration  equipment,  yielded  results  within  0.1  °C  of  the  reference 
temperature. 
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Table  2.  Platinum  Wire  Calibration  Data  (d  =  127  )xm) 
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Temperature  (C) 


55 


Temperature  vs.  Resistance 
(50.8  pm  diameter) 


Figure  15.  Platinum  Wire  Calibration  Curve  (d  =  50.8  pm) 


46 


48 


APPENDIX  B:  SAMPLE  CALCULATION 


The  following  sample  calculations  use  data  collected  during  an  experimental  run 
at  f  ~  1213  Hz  (L  =  64  cm),  pressure  ratio  ~  1.0%  and  AT  ~  24°C  for  the  50.8  pm 
diameter  platimun  wire.  The  measure  values  and  the  properties  for  air  (at  300K)  are 


=25.3‘C 

/  =  U\%Hz 

=  653.18mF 

Vpf,  =  14.03 10/M  F 
F„  =  0i33F 


k  =  0.026'%,jf 

v  =  li9xl0‘*"*% 
P„  =  101.35itPa 


where  Vn,  is  defined  as  the  measured  pressure  transducer  output  voltage  (RMS)  and  Pm  is 
the  mean  ambient  air  pressure. 

A.  NUSSELT  NUMBER 

The  Nusselt  number  equation  is 


where 


Substituting, 


Power 

A.(t,  -  T.) 


Nu  = 


Power  ( 


(B.l) 

(B.2) 


(B.3) 
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The  wire  surface  temperature,  Ts,  is  calculated  using  the  equation  found  from  the  wire 
calibration  as  described  in  Appendix  A. 


and 


7;  =  59.8156/2,, -224.2799 


V 


(653.1 8m  F)(98.3mQ) 


(14.03 10m  F) 


Therefore, 


=  4576.1 2mn 


r  =49.4'C 


Power  is  found  from 


Vj  (653.18mF)f^X00mFF 

K  (4576.12mn)(l%oo;;,Q) 


(B.4) 

(B.5) 


(B.6) 


Power  =  0.0931F 


Substituting  into  the  Nusselt  number  equation. 


_ (0.0931F)2(254xl0~^m) _ 

(1.20x10-' m'X49.4  -  25.3)(^)(o.026%j^) 


where 

4  ^ndl  =  2;r(2i4xl0-'m)(7i3xl0'^ m)  (B.7) 

4  =  1.20x1 0-'m^ 


50 


Finally, 


Nu  =  0.628 

B.  REYNOLDS  NUMBER 

The  Reynolds  number  equation  is 


„  USla)  ^  , 

Re^  =  =  Ish^ 


where 


Uo  =  - 

y 


^  p\ 

f± 


c  =  ^r^(r„  + 273.15) 


and 


P_ 

"  0i22 


(B.8) 

(B.9) 

(B.IO) 

(B.ll) 


Making  the  above  substitutions  into  the  Reynolds  number  equation,  we  get 


Re^  = 


0.522  vR 

ffl 


(B.12) 


Re^  — 


2(2.54x1 0~^/w)(0.533F) 


(o.522  1.59x1 0-'  01.35yfc?a) 


1 


287 


m 


1.4 


^V(25.3  + 273.1 5)(^) 


Finally, 
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APPENDIX  C:  UNCERTAINTY  ANALYSIS 

A.  PLATINUM  WIRE  RESISTANCE 

The  equation  for  platinum  wire  resistance  is 


where 


7  = 


therefore, 


The  uncertainty  becomes 


R 


R  \ 

W  W  jr 


'^R,wire 


*V,wire 


■  ^R,PR 


PR 


\^PR 


V 


H’,PR 


(C.l) 

(C.2) 

(C.3) 

(C.4) 


^R,wire 


+ 


K 


(iMi 

I  f',/ 


Y 

UypR 

/ 


(C.5) 


u, 


V,wire 


V  K 


+ 


R 


+ 


-U, 


V,PR 


V  ypR  J 


(C.6) 


7?;,^  =  0.0983/2 
=  0.026% 

=  measured  quantity 
^p.wirc  =  reading  +  #  counts) 


Resistance,  Precision  Resistor 
Uncertainty  in  Resistance,  Precision  Resistor 
Voltage,  Platinum  Wire 
Uncertainty  in  Voltage,  Platinum  Wire 
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Vf,^  =  measure  quantity  Voltage,  Precision  Resistor 

Uy  PH  =(%  reading  +  #  counts)  Uncertainty  in  Voltage,  Precision  Resistor 


The  measured  quantity  uncertainty  equations  vary  depending  upon  the  Hewlett 
Packard  HP  3478A  digital  multimeter  range  and  resolution  for  each  specific  data  sample. 
The  following  is  an  example  of  a  data  run  for  f  ~  1213  Hz  (L  =  64  cm),  pressure  ratio  ~ 
1.0  %  and  AT  ~  24°C  using  the  50.8  pm  diameter  wire. 


Rpp  =  0.0983/2 


F=  653.1 8wF 


Up  pp  =  0.026% 


Vpp  =14.0310/wF 


r  0.006%' 


*V,wire 


X 


100  J 


653.1 8wF  +  Icountl 


^100jdO-‘’FVlOO0wF 


^  count  J 


V  F 


=  0.13919/wF 


^V,PR  ~ 


("0.035%^ 

r,  flOOxlO^’Fl 

14.030mF  +  A0counts\ - 

100  J 

V  count  y 

lOOOmF^ 

F  J 


Uy  pp  =  0.0089 1/wF 


The  platinum  wire  resistance  uncertainty  is 

Up^_  Ir0.13919y  r  0.026 Y  ^  T -0.0089  iV 
”vl  653.18  J  100  J  14.030  J 
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-^  =  7.186x10-'  or  0.07186% 

R.. 


The  calculated  platinum  wire  resistance  is 


R^  =  ASieXlmQ 


therefore. 


^RM,re  =  X29mn 


B.  NUSSELT  NUMBER 

The  Nusselt  number  equation  is 


(C.7) 


where 


therefore, 


Power 


w  / 

/R.. 


RJ  (d 


dl7c{T,-TjVk)  RJk7r(T,-TJ 


From  the  calibration  data  the  surface  temperature  of  the  wire  can  be  written  as 


T.  =  AR^,  +  B 


(C.IO) 


where  A  and  B  represent  the  slope  and  y-intercept,  respectively,  from  the  linear 
regression  model  used  for  the  wire  calibration.  Substituting,  the  Nusselt  number 
becomes, 
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m  = - ^ - —  (til) 

lkxRJAR.+B-TJ 

and  the  uncertainty  becomes, 


The  following  continues  from  the  previous  example  where  f  ~  1213  Hz,  pressure 
ratio  ~  1 .0  %  and  AT  ~  24°C. 

^  =  59.8156°C/Q 
5  = -224.2799' C 
T^=25.TC 
Uj  =  0.5'  C 
=  4576.1 2/W/2 

UR..,rc  = 

The  Nusselt  number  uncertainty  is 

^  =  ^(4.26x1 0-')^  +(-8.90x10-')^  +(2.08x10-^)^  +(3.00x10’')^ 
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=653.1 8wF 
=0.13919mF 
/  =  0.0753m 
u,  =  3% 


—  =  3.75x10-^  or  3.75% 
Nu 


The  calculated  Nusselt  number  is 


Nu  =  0.628 


therefore, 

=  0.0236 


C.  REYNOLDS  NUMBER 

The  Reynolds  number  equation  is 


where 


Re^  = 


{2a)U, 


=  2- 


a 


£  4 
y/pj 


(C.13) 


c  =  ylrR{T^+ 27315) 


(C.14) 


and 


(C.15) 


V2 


=  microphone  output  (RMS) 


Substituting, 


Re^  = 


R  \ 


1 


0.522  P„ 


(2a)[^j^(7; +273.15) 


(C.16) 


Let 
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(C.17) 


Fo'  =  ^  Volts ,  and  =  (7^  +  273.1  S)K 


The  uncertainty  becomes 


«Re  = 


r^Re. 


cti  “ 


+ 


^7Re^  I  ^ 


V  “""7 


(C.18) 


«Re 

2 

+ 

fu  ■) 

2 

+ 

'‘r.  ] 

Re,  V 

\  aJ 

[Vo) 

[2Tj 

(C.19) 


The  following  continues  from  the  previous  example  where  f  ~  1213  Hz,  pressure 
ratio  ~  1 .0  %  and  AT  ~  24°C. 


a  =  2.54x10 
u  =  5% 


7;'=  (25.3 +  273.15)/: 


M  .  =  0.5/: 

^ac 


F,  =  0.533F 


,  ro.06%^  fO.03%' 

u„.  =2.3xlO“V+  -^tt-|0.533F  +  1— t:-|1F 


100 


100 


u...  =  2.920x1 0'^ 


The  Reynolds  number  uncertainty  is 


58 


u 


Re 


Re^ 


f  2.920x10-'^ 

2 

+  1 

1  0.533  ; 

1 2(298.55)  J 

=  0.0503  or  5.03% 
Re, 


The  calculated  Reynolds  number  is 


Re,  =  7.96 


therefore, 
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APPENDIX  D:  EXPERIMENTAL  DATA 


This  appendix  includes  a  representative  sample  of  experimentally  measured 
quantities,  calculated  parameters  and  uncertainties  for  the  two  platinum  wire  diameters 
used.  Due  to  space  constraints,  some  constants  used  in  the  spreadsheets  for  calculations 
have  been  omitted.  The  50.8  pm  diameter  wire  data  is  reported  first  followed  by  the  127 
pm  diameter  wire. 
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iai8 


63 


6.02 


67 


536. 


73 


5.02 


75 


77 


46.55 


,30E'9 


79 


08 


81 


38 


83 


85 


29.8  2.08  0,001  8.90  28.0  1.67  1.07  3.81  E-7  101.3 
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lE-4 


30.83  8.23  440  360  1049  0.901  168  501 
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2.02 


% 


06 


91 


663 


>56.8  562 
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